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We give a theoretical description of avalanche- like dynamics of magnetic flux in the critical state 
of "hard" type-II superconductors using a model of a one-dimensional multijunction SQUID that 
well reproduces the main magnetic properties of these objects. We show that the system under 
consideration demonstrates the self-organized criticality. The avalanches of vortices manifest them- 
selves as jumps of the total magnetic flux in the sample. The sizes of these jumps have a power-law 
distribution. Our results are in qualitative agreement with experiments. 
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I. INTRODUCTION 

About three decades ago Charles P. Bean proposed 
a model describing a critical state in "hard" type-II 
superconductors." In the frame of this model, the critical 
state, arising in a superconductor as a result of inter- 
play of magnetic and pinning forces, is a non-equilibrium 
state, in which the density of Abrikosov vortices varies 
linearly with the distance from the sample surface. Later 
Pierre G. de Gennes considered an interesting analogy 
between the critical behaviour of magnetic flux in su- 
perconductors and evolution of a sand pile.^ De Gennes 
compared the magnetic flux dynamics in the critical state 
with the process that is observed in an over-critical sand 
pile where sand slides from the top of the pile forming 
avalanches. According to this brilliant concept, when the 
superconductor becomes over-critical, the vortices begin 
to move in form of avalanches. Each avalanche includes 
a simultaneous motion of a large number of vortices. 

Although this picture of the magnetic dynamics in su- 
perconductors was formulated in 60-th, the experimental 
technique, allowing observation of vortex avalanches and 
studying their size distribution, was developed only in the 
last decade (see, for example, Ref. 0) . These experiments 
demonstrated that in certain conditions slow changes of 
an applied magnetic field result in an avalanche-like re- 
distribution of vortices. It turns out that the sizes of 
avalanches have a power-law distribution^*^ Such a dis- 
tribution is, as is well known, a characteristic feature of 
systems with a self-organized criticality (SOC)ii 

The concept of SOC was formulated in 1987 in order to 
explain a behavior of giant dissipative dynamical systems 
consisting of a large number of interacting elementsJ- Ac- 
cording to the main principles of this concept, such sys- 
tems are able to accumulate small external perturbations 
and evolve into a self-reproducing critical state. This 
state is an ensemble of various metastable states. The 
critical system migrates from one metastable state to the 
other by means of avalanches. The avalanches are initi- 
ated by small external perturbations and may have vari- 
ous sizes. The distribution of avalanche sizes is a power- 
law function with an exponent close to 1. The paradigm 



of SOC is a sandpile. This is why the corresponding 
mathematical model is called the sandpile modeliSi 

Recently, considering the similarity between the criti- 
cal magnetic dynamics in superconductors and SOC, we 
showed that employing of dynamical equations provides 
a promising way to correlate these two phenomenaAiSiii 
In these works, we used the fact that the main features 
in the behavior of a "hard" type-II superconductor may 
well be reproduced by considering a so called discrete su- 
perconductor, i.e., a lattice of Josephson junctions (mul- 
tijunction SQUID) iSiiS An additional advantage of this 
approach is that the corresponding equations are not dif- 
ficult for the analysis. 

It is known that magnetic properties of a lattice of 
Josephson junctions may be characterized by a param- 
eter V ~ jcCL'^/^o {jc is the critical current density of 
the junction, $o is the magnetic flux quantum, and a is 
the interjunction distance) liSiH In the case oi V <C 1, 
the system can be considered as a single sandwich-like 
Josephson junction without pinning.^'* If the condition 
y ^ 1 is not satisfied, a multijunction SQUID, due to 
its discreteness, is able to pin magnetic flux. In the case 
of y 3> 1, we have a system with strong pinning. 

It was also shown theoretically and by computer sim- 
ulations that the critical state of two-dimensional (2D) 
and one-dimensional (ID) multijunction SQUID's with 
y ^ 1 is self-organizedjSiiMi Avalanches in such sys- 
tems manifest themselves as pulses of voltage across the 
j unctions It was demonstrated that the total voltage 
integrated over the time of the avalanche is an analog of 
the avalanche size in the sandpile model. ^ For some way 
of perturbation, this quantity has a power-law distribu- 
tion. It was also proven that for V ^ 1 the dynamical 
equations for the superconductor are equivalent to algo- 
rithms describing the sandpile model. 

The quantities usually measured in the sandpile exper- 
iment are the mass of the pile and its variations,— The 
direct analog for the multijunction SQUID is the total 
magnetic flux $ in the sample and its fluctuations. This 
characteristic was studied experimentally, for instance, in 
Refs. 0and0. This is why, in this work, we concentrate 
on calculation of 

The aim of this work is a theoretical investigation and a 
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computer simulation of the magnetic flux dynamics and 
the avalanche statistics in the critical state of a multi- 
junction SQUID. We use a model of a ID multijunction 
SQUID with intrinsic spatial randomness, which was in- 
troduced in Rcf. JJ, This model takes into account that 
real superconductors are disordered systems. A spatial 
randomness, introduced in our model, is an equivalent of 
the existence of randomly distributed pinning centers^. 
We consider here not only the case of U 3> 1, studied 
earlier jSiiSiii but also F < 1. We show that for all values 
of V and for a fixed degree of disorder, the critical state 
of the system may be represented as a set of metastable 
states transforming to each other due to avalanches. An 
avalanche in such a system results in penetration of mag- 
netic flux into superconducting sample and its redistri- 
bution between the system cells. We demonstrate that 
for all considered values of V the probability density of 
magnetic flux jumps has a power-law distribution. Thus, 
using the simplest model of a superconductor, we obtain 
the results that are in qualitative agreements with exper- 
iments. 

The paper is organized as follows. The second section 
is devoted to description of our model of ID multijunc- 
tion SQUID with intrinsic spatial randomness. We ana- 
lyze the dynamical equations and demonstrate that dis- 
order of the system is a key factor for realization of the 
complex critical state. In the third section, we discuss 
the structure of the critical state for various values of the 
SQUID-parameter V. Single avalanches are described in 
the forth section. In the fifth section, we analyze the 
statistics of magnetic flux jumps and discuss the possi- 
bility of the realization of the self-organized critical state 
in the system under consideration. The results of calcu- 
lations are compared with experiments. In Conclusions 
we summarize the main results of the work. 



II. BASIC EQUATIONS 

The one-dimensional multijunction SQUID, which we 
consider here, may be represented as two infinitely long 
in the y-direction superconducting layers connected by 
Josephson junctions, as is shown in Fig. ^ All junctions 
have the same length / along the x-axis. The distance 
between the i-th and the i + 1-th junctions we denote as 
a random variable bi. The system is placed in a slowly 
increasing external magnetic field iJcxt , aligned along the 
y-axis. Using the resistive model of a Josephson junction 
and neglecting by thermal fluctuations, we can write the 
current density ji as: 



ji = jc sin ipi + 



$0 dip^ 
2ttp dt ' 



(1) 



where jc is the critical current density, ipi is the gauge- 
invariant phase difference across the i-th junction, p is the 
junction resistance per unit area. The current density ji 
is connected with the magnetic fleld Hi in the neighbor- 
ing cells by the following expression (we numerate the 
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FIG. 1: The x — z section of ID multijunction SQUID 



cell by the nearest left junction): 
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where = HiSi is the magnetic flux in the i-th cell. Si — 
2Xbi is the cell area, A is the magnetic field penetration 
depths 

Taking into account that 
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[iPi+i{t) ~ <fii(t)] 



(3) 



we obtain the system of differential equations for the 
gauge- invariant phase differences: 



V sin ipi + T 



dt 



[Jifi+i 
V sin If I + 

V sin Lpj^ -\- T 



{Ji + Ji^i)tfi + Ji-ifi-i] ; i ^ 1, iV; 
[Jiip2 - Jwi] - 27r/ioxt; 
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where a — {bi) is the average interjunction distance. 

In order to take into account a spatial disorder of real 
physical systems, we consider the distances between the 
junctions as random. As may be seen from Eqs. Q), in 
our model this randomness is equivalent to a scatter of 
the coefficients Ji. In the following, we assume that the 
values of these coefficients are distributed uniformly in 
the interval [1, 1 -I- A J]. 

As was stated above, there must be a large number 
of metastable states in the system in order to expect 
a self-organized behavior. In our case, such a situa- 
tion can be realized only if AJ^ ^ 0. This is illus- 
trated in Fig. [3 which shows the dimensionless current 
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FIG. 2: The dependence of the dimensionless current z[k] for 
fc = 3 and V = 0.6 on the number of avalanches n for A J — 
(upper panel) and A J = 0.1 (lower panel) scattering of the 
coefficients Ji. We see that introducing of disorder results in 
a crucial increase of the number of metastable states. 



z[k] ~ y/27rsin(^fe in one of the junctions at the final 
moment of rt-th avalanche for two cases: (i) for a regular 
system with AJ; = and (ii) for a disordered system 
with AJ = 0.1. We see that, while the regular system 
behaves itself periodically and has only a limited number 
of possible values of the junction current, the situation 
for AJ ^ is completely different. In the latter case, 
the number of possible values of Zk is rather large. This 
result does not depend, of course, on the number of the 
chosen junction. 



III. STRUCTURE OF THE CRITICAL STATE 
OF ID MULTIJUNCTION SQUID 

In this section, we present the result of computer simu- 
lation of the critical state for the system described above. 
Four different values of the SQUID-parameter are con- 
sidered. We take F = 40 in order to model the case of 
y > 1 and y = 0.3, V = 0.6, and V = 1.2 to simulate 
the situation with < 1. 

We use an Euler integration scheme for Eqs. Q with 
Ai = 0.01, for y > 1 and with Ai = 0.1 for V <\. The 
non-stationary method of perturbation, which is com- 
monly used for simulations of self-organized systems, is 
employed. It means that the external magnetic field /loxt 
varies only when all relaxation processes in the system 
are already finished. This way of changing of the exter- 
nal field is close to that in experiments of Aegerter et 
aU 

Before starting, we fix a set of random coefficients 



Ji, which remain unchanged during the simulation pro- 
cess. Starting from the state with ipi — 0, we perturb 
the system by increasing the external field from /icxt to 
/loxt + A/ioxt- In our simulations, we use A/ioxt = 1 for 
y = 40 and A/icxt = 0.1 for small and transient values of 
V . Then we allow the system to relax and, as was men- 
tioned above, we assume that the value of /icxt is constant 
during the relaxation process. We take that the system 
has already reached its metastable state if dipi/dt < 10" 
for all i. When the dynamics stops we perturb the system 
again and so on. 

First, we analyze the distribution of magnetic field in- 
side our SQUID for various values of V. The magnitude 
of the dimensionless magnetic field in the i-th cell, mea- 
sured at the end of n-th avalanche, may be calculated as: 



2Aa$ 
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(5) 



Note that for all values of V and A J the system demon- 
strates irreversible magnetic behavior. Remanent mag- 
netization becomes zero only for negligible value of the 
parameter V {V — 0.06) and for A J^ — 0. 

Fig. 13 shows the profiles of magnetic field inside the 
SQUID for one of metastable states for two different val- 
ues of V and A J = 0.1. We see that for V = 40, when 
every cell acts as pinning center, the profile is similar to 
the result of Bean's model. In the case of ^ = 0.3, the 
pinning centers are represented by groups of cells and a 
number of peaks may be seen. We note that the peak 
amplitudes are different for different metastable states. 

Now we consider the dynamics of the total magnetic 
fiux in the sample and its dependence on the external 
magnetic field /icxt- The total magnetic flux may be cal- 
culated as 



Here t/?^"'* denotes the value of the phase for the final 
moment of n-th avalanche. The variation of the total 
magnetic flux due to the n avalanche 

A$(") = $(") - $("-1). (7) 

Figures 01a) and^Jb) illustrate the evolution of the total 
magnetic fiux with increasing magnetic field for two sys- 
tems with different values of V. As was noted earlier, the 
external magnetic field is assumed to be constant during 
the avalanche. The values, plotted in Fig. 01 correspond 
to the moments when the relaxation is already finished 
but a new step of the external field is not yet applied. 
Fig. ^Jc) and^Jd) show the corresponding dependencies 
of A<I>. As may be seen, there are avalanches of different 
sizes. We also see that the behavior remains qualitatively 
the same for different values of V. 

In Fig. 13 we show experimental results of Ref. 5. In 
this work, the dynamics of the magnetic flux in the crit- 
ical state of a thin YBa2Cu307_3. fllm was studied. The 
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FIG. 3: The magnetic field profiles h]"' inside the SQUID for 
V — 40 (upper panel) and V = 0.3 (lower panel), A'' = 257. 




FIG. 4: The total magnetic flux in the SQUID (a-b) and the 
magnitude of the flux jumps (c-d) as functions of the external 
magnetic field. The SQUID size iV = 65 and AJ = 0.1. 
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FIG. 5: A$ as a function of time from Ref. |^ The data pre- 
sented in this figure come from series of nine experiments and 
different experimental runs are separated by vertical arrows. 

external magnetic field was changed in step- wise manner. 
After each field step of 0.5 Oe the sample was allowed to 
relax for 10 seconds before the magnetic flux was mea- 
sured. These experimental conditions arc similar to rules 
that we use in our computer simulation. As may be seen 
in Fig. [SI similar to that in computer simulation, ex- 
perimentally measured magnitudes of the magnetic flux 
jumps have various sizes. 

The magnetic flux jumps was also observed in an arti- 
ficial 2D lattice of Josephson junctions.^ In these experi- 
ments, the external magnetic field was changed continu- 
ously but slowly. This allows to neglect by the change of 
the magnetic field during the relaxation process. Fig. El 
shows the hysteresis loops of the total magnetic moment 
of the lattice from Ref. |^ It may be seen that the 
avalanches, which are represented by jumps of the mag- 
netic moment, have different sizes and include hundreds 
of magnetic flux quanta. The upper trace shows several 
superposed hysteresis loops for the same experimental 
conditions. Due to this superposition, we can see the 
randomness in jumps. The random jumps with differ- 
ent sizes may be considered as a manifestation of SOC. 
As may be seen comparing Figs. 0] and the results of 
our simulation are in qualitative agreement with experi- 
ments. 



IV. THE AVALANCHE STRUCTURE IN THE 
CRITICAL STATE 

As was demonstrated above, an increase of the external 
magnetic field can launch an avalanche. In this way the 
system migrates form one metastable state to the other. 
An avalanche in our multijunction SQUID is a simulta- 
neous penetration of a considerable number of vortices 
in the sample and a redistribution of the corresponding 
magnetic fiux inside the system. As a result of such a 
process, the system reaches a next metastable state, the 
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FIG. 6: The histeresis loops of the total magnetic moment in 
the 2D lattice of Josephson junctions from Ref. 0. 



total magnetic flux increases and the values of the mag- 
netic flux in some cells change. In this section, we con- 
sider the process of the avalanche development. 

At every moment i„ during the n-th avalanche we 
calculate the magnetic flux ^'f'\tn) using the expres- 
sion (jSJ. We also calculate the difference A$j-"''(i„) — 

(i„) — $i"''(ion)] where ton is the initial moment of 
the avalanche. We emphasize that all the results in this 
section are related to a single avalanche and the external 
magnetic field is assumed to be constant. 

Fig. Uldemonstrates the process of penetration of mag- 
netic flux inside the SQUID and its redistribution be- 
tween the cells during the avalanche. The results pre- 
sented in Fig. correspond to a regular SQUID with 
A J = 0. The flgure shows the magnetic flux distribu- 
tions corresponding to several different moments in time 
during the avalanche. We see that at the beginning of the 
avalanche the magnetic flux penetrates into the bound- 
ary cells and only with time the penetration reaches the 
central part of the SQUID. The analogous results for a 
disordered system are presented in Fig. Eb- In this case, 
the magnetic flux demonstrates the same dynamics as 
for a regular system. However, due to the randomness 
the magnetic flux proflle is much less regular than in the 
previous case. 

li V ^ 1, each flux quantum may be considered as 
distributed between several neighboring cells. In some 
cases, motion of such extended vortices results in negative 
values of A$, as may be seen in Fig.|Sl This is true for 
both regular and disordered SQUID's. 

V. THE AVALANCHE STATISTICS IN THE 
CRITICAL STATE AND THE SELF-ORGANIZED 
CRITICALITY 

In experimental works on the critical state of super- 
conductors, the authors often points out on the similar- 
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FIG. 7: The distribution of the magnetic flux in different 
moments during a single avalanche. V = 0.6, A'^ = 40. Time 
is measured in arbitrary units, (a) AJ — 0; (b) AJ = 0.1. 



ity of the system behavior and the phenomenon of self- 
organized criticality. Besides the avalanche- like dynam- 
ics, it is a power-law distribution of magnetic flux jumps. 
For instance. Fig. |^ shows the distribution of avalanche 
sizes from Ref. A. In this work, the avalanche-like dynam- 
ics of magnetic flux in a NbTi tube was studied. Three 
distribution functions correspond to three flxed values of 
the external magnetic held. The external magnetic field 
varies in the interval of 30 Oe centered at one of the fixed 
values with the rate of 5 Oe/s. 

Fig. ^1 demonstrates the results of our calculations for 
the probability densities of jumps of the total magnetic 
flux for multijunction SQUID for three different values 
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FIG. 8: The magnetic flux differences A^\n){t„)/$o during 
a single avalanche. V = 0.6, A'^ = 65 and A J = 0.1. 
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FIG. 10: Probability densities for the magnetic flux jumps 
p(A$/$o) for AJ = 0.5. (a) V = 40, (b) V = 1.2, (c) V = 
0.6. The straight lines represent a power law: (A$/$o)" with 
a — —1.01, —1.26, —1.59 for (a), (b) and (c), respectively. 



ical state, in which the total magnetic flux plays a role 
of the main characteristic of the system. 
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FIG. 9: Distribution functions of the magnetic flux jumps 
obtained in Ref. |^ 



VI. CONCLUSIONS 



of the parameter V. As was shown earlier, "'^-'^^ the case 
of F = 40 for this degree of disorder demonstrates the 
self-organized behavior. As may be seen in Fig.^| there 
are rather extended intervals of power-law distributions 
for all considered values of V. 

At the same time, from the point of view of classi- 
cal interpretation of SOC the magnitude of flux jumps is 
not a direct analog of the avalanche size.'' However, just 
as the classical self-organized critical state, the critical 
state in real superconductors, as well as in our model, 
is self-reproducing and consists of a large number of 
metastable states that transform to each other by means 
of avalanches. Thus, we can conclude that here we are 
dealing with a more general type of a self-organized crit- 



Based on dynamical equations describing the simplest 
model of a discrete superconductor (ID disordered mul- 
tijunction SQUID), we present a theoretical description 
of the avalanche-like dynamics of the magnetic flux in 
"hard" type-II superconductors. Different values of the 
SQUID-parameter V ^ jca^/^o are considered. For all 
values of V, including F < 1, the critical state in the 
multijunction SQUID can be considered as a general- 
ized type of a self-organized critical state. In contrast 
to the classical definition of SOC,'' the main characteris- 
tic of generalized critical state is a size of magnetic fiux 
jumps. This quantity demonstrates a power-law distri- 
bution if some degree of disorder is introduced into the 
system. Our results are in qualitative agreement with 
experiments. 
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